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Abstract 

The transition from the low Q 2 region in deep inelastic scattering where charm produc- 
tion is described by photon-gluon fusion to a region where the charm structure function 
F% is largely generated by a charm quark density should be a smooth passage. The all 
orders prescription for matching between the two schemes does not uniquely determine 
all the relevant coefficient functions (CF's) or guarantee the correct threshold behaviour 
at fixed order in a s . The constraints of matching not only the value but also the evolution 
of F£ order by order allows all CF's to be determined at any order, each of which reduces 
to the appropriate massless MS expression as Q 2 — > oo, and leads to a smooth transition 
in the threshold region. This procedure gives an excellent description of the F£ data for 
all Q 2 . 
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1 Introduction 



Recent measurements of charm production at HERA |], § emphasise the importance of hav- 
ing a consistent theoretical framework for heavy flavour production in deep inelastic scattering 
(DIS). The charm structure function F 2 C could be 20% or more of the total F2 and a modern 
global analysis of structure functions must necessarily include a satisfactory description of F£. 
Recently we have developed a general order-by-order procedure which solves the problem of 
how to calculate F| in the region close to Q 2 = m 2 where m c is the charm quark mass. Below 
Q 2 = m 2 the conventional description in terms of order-by-order in a s boson-gluon fusion is 
perfectly adequate, but for Q 2 > m 2 this description becomes increasingly unreliable due to 
the presence of potentially large logarithms at all orders in a s which should be resummed. By 
changing to the alternative description where the charm quark is treated as a parton, at the 
point Q 2 = m 2 , this resummation is automatically performed and, at the same time, a com- 
plete set of parton densities needed to calculate other processes involving nucleons is obtained. 
While this idea of transferring from one description to another is not new, until now there has 
not been a satisfactory procedure for preserving threshold behaviour across the transition point 
whilst systematically presenting the calculation to any well-defined order. 

In this paper we present our results at leading order (LO) and next-to-leading order (NLO) 
to demonstrate the key features of our method, the practical advantages it offers and how it 
improves upon previous attempts to treat heavy flavours. We demonstrate that the require- 
ment of all orders equivalence between the two descriptions does not actually provide unique 
expressions for every relevant CF. The order-by-order in a s matching of the structure function 
itself between the two descriptions at the transition point does not remove the ambiguity. Extra 
information is required to do this and, in our case, this is provided by the further requirement 
that the order-by-order evolution of the structure function should also be continuous across 
the transition point. This added constraint then allows each CF to be determined with the 
practical benefit that a very smooth transition is ensured, the correct threshold behaviour being 
maintained across the transition point and above. Here we only sketch how our results are de- 
rived, the emphasis here being on how the procedure works in practice; the detailed theoretical 
justification is presented in full in ||. 



Consider the case where charm is produced by the neutral current. At low Q 2 , the number of 
active flavours, rif, is 3 and the charm cross section is generated by photon-gluon fusion (PGF) 
0, |5|. This corresponds to the so-called fixed flavour number scheme (FFNS) and provides an 
acceptable description provided Q 2 is not large and one is not interested in the concept of a 
charm quark density. 

In the FFNS at order 0(a s ) the charm structure function is given by 
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The FFN, ZM-VFN and VFN schemes 



F^x,Q 2 ,m 2 c ) 
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C^ FF (Q 2 /m 2 c ) ® ^ /=3 (/i 2 ), 



(1) 
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where the CF which is convoluted with LO evolved gluon density g nf =3 is 

1 + v 



C (l) FF { 



z, e 



(Pq g (z) + 4ez(l - 3z) - 8eV) In 



+ (8z(l - z) - 1 - 4ez(l - z))v 



l-v 

2 /| m 2 



6{W 2 - 4m 2 c ) (2) 



where e = m 2 /Q 2 , W 2 = Q 2 (l/z — 1), the gluon quark centre of mass energy, v is the velocity of 
the charm quark or antiquark in the photon-gluon centre-of-mass frame, defined by v 2 = 1 — 
Am 2 /W 2 , and Pg g (z) = z 2 + (l — z) 2 , the LO quark-gluon splitting function. These f-dependent 
terms ensure that the coefficient function tends to zero smoothly as W 2 = Am 2 is approached 
from below, and hence the structure function has a smooth threshold in W 2 = Q 2 (l/x — 1). 

In the limit Q 2 — > oo the gluon CF reduces to 

*F( Z , e ) P° g (z) \n(e) + P» g {z) In(^r^) + 8z(l - z) - 1 (3) 

from which we can see the potentially large logarithm lne = ln(Q 2 /m 2 ). The natural way to 
sum this contribution together with higher powers of such logs is to describe charm production 
through a charm density which evolves according to the standard (light quark) renormalisation 
group equations. 

This leads us onto an alternative description which is simple (and therefore practical) and 
which is commonly used in global analyses of DIS data || [7| where the aim is to extract parton 
densities for both light and heavy flavours. In this scheme, which we label the zero mass 
variable flavour scheme (ZM-VFNS), a value of [i 2 (~ m 2 ) is chosen below which m c is taken 
to be infinite and above which m c is taken to be zero. While this is intended to be a rough 
description only it has proved successful in describing the then existing (pre-HERA) charm 
structure function data well above the threshold region with fx 2 w 3GeV 2 . However it is bound 
to be unrealistic at low Q 2 because of the abrupt variation in Q 2 - in contrast to the proper 
smooth threshold behaviour in W 2 . As Q 2 increases the ZM-VFNS becomes more reliable, the 
0(a s ), i.e. the NLO expression in this case, being given by 

F^x,Q 2 ,m 2 c ) = C:f=\Q 2 /^ 2 ) ® c+t/Am^ 2 ) + C n g f=\Q 2 /^ 2 ) ® g nf= ^ 2 ,m 2 J^ 2 ) (4) 

where the charm density c + = c + c and the CF's are expanded as 

CT'= 4 (z,l) = C?\z) + ^-C^{z) 

c;f=\z,i) = %cp m (z) (5) 

where 

C<°\z) = z 5(1 - z) and Cf )m = P° qg {z) ln(— ^) + 8z(l - z) - 1. (6) 

Our aim is to develop an approach which extrapolates smoothly from the FFNS at low Q 2 
to the ZM-VFNS at high Q 2 , maintaining the correct ordering in both schemes. First we note 
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that in the FFNS eq(@) is valid up to corrections of 0(A 2 /m 2 ) while in the ZM-VFNS eq(g) is 
valid only up to corrections of 0(m 2 /fj, 2 ). In order to improve the accuracy of eq(|]) we need 
to examine the connection between the parton densities in the two schemes. Eq(|]) actually 
has an apparent extra degree of freedom over eq(^), namely the charm density evaluated at 
its starting scale. This is of course not really true in practice and the connection between the 
parton densities for 3 and 4 flavours takes the form 

c + (z^\m 2 Jii 2 ) = A c \ii 2 /rr%) ® g nf =^ 2 ) 
g nf=4 (z,fx 2 ,m 2 c /fx 2 ) = A 99 {^/m 2 c ) ® g nf=3 (fx 2 ) (7) 

at leading order, where the elements A ba which contain ln(/i 2 /m 2 ) terms, are, in general, part 
of a full 5x4 matrix which also connects the light quark flavours. 

For Q 2 3> m 2 , the equivalence of the FFNS and the ZM-VFNS at all orders then lead to 
the connections between the CF's in the two schemes up to 0(m 2 /[i 2 ) ||, in particular up to 
0(a 2 ) 

C FF (z,Q 2 /^ 2 ,Q 2 /ml) = C n J=\Q 2 lf) ® A"(f/ml) 

+ C"/= 4 (Q 2 //i 2 ) ® A^/m 2 c ) + 0{m 2 J^ 2 ). (8) 



The details of the connection are fully worked out in [[|. To improve the accuracy of 
where the uncertainty is reduced to 0(A 2 /m 2 .), requires defining 'corrected' CF's, C^ F , in 
another nf = 4 scheme - the variable flavour number scheme (VFNS) - where one can write 

+ Cj F (Q 2 /fi 2 ,m 2 c /fi 2 ) ® ^ /=4 (/i 2 ,m 2 //i 2 ) + (9(A 2 /m 2 ), (9) 
where the corrected CF's are related to the FFNS CF's by 



(10) 



the new n/=4 CF's now being exact at all values of Q 2 . 

In the spirit of ref|| the procedure we adopt is to use the FFNS for Q 2 < m 2 and to switch 
to the VFNS for Q 2 > m 2 . The precise choice of the transition point is undetermined, though 
is naturally of order m 2 . However, there are matching conditions between the partons in the 
two schemes in order to ensure consistency with the usual asymptotic expressions, e.g. eq(7), 
and taking yU 2 = m 2 removes complications arising from ln(/x 2 /m 2 ) terms in these matching 
conditions. To make practical headway in the VFNS one must solve eq(^) for the C^ F . 
However one can see that the all orders matching of F% in the two schemes, from which eq (|T0"|) 
arose, is not sufficient since, for example, at low orders the single quantity C FF is expressed in 
terms of the two quantities C^ F and Cj F . We stress that any choice satisfying eq(8) is "correct" 
in the sense that it leads to the same all orders expression. Nevertheless, each choice leads to 
a different expression if one uses the usual rules of combining coefficient functions and parton 
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distributions of a given order to obtain a fixed order in a s expression for the structure function. 
In principle a completely consistent theoretical approach would lead to a completely seamless 
transition between the two schemes at every order in perturbation theory. In practice this seems 
to be, at the least, incredibly difficult and is probably not even possible in any unique manner. 
However, in order to obtain as well-defined a theoretical procedure as possible we impose not 
only continuity of the structure function but also demand, in addition, order-by-order matching 
of the evolution of F| in the two schemes, as we shall see below. 

3 The practical solution 

The explicit form of eq(0) for /i 2 = Q 2 at 0(a s ) is 



i2\ _ „ /_ rSi 



® 9n f =3 



g nf =ii.z } Q ) = g nf=3 (z,Q ) - — In — g n 3 . (11) 

Inserting the implied expressions for the matrix elements A C9 (z, /i 2 /m 2 ) and A" \z , /i 2 / m 2 ) into 
eq(8) gives the familiar relation 

Cf FF (z, Q 2 /m 2 c ) = C« VF (z, Q 2 /m 2 c ) + Cf VF (Q 2 /m 2 c ) ® P° ln(g) (12) 

connecting the gluonic CF's in the FFNS and VFNS. In previous implementations of a VFNS 
|10| , eq([T2|) served as the definition for VF in terms of the PGF CF (eq@) with an assumed 
form of Cf^ VF given by 

VF (z, Q 2 /M 2 ) = z 6(x - Z )(l + ^, x = (l + ^) _1 (13) 

where the delta-function describes the tree-level diagram for a massive quark scattering from 
a photon and the modified argument of the delta-function follows from demanding that the 
massive quark is on-shell. We note that the same definition of the zeroth order coefficient 



function is adopted in [II | though there are differences between this and [IT]], notably a mass 
dependent evolution. 

We do not believe that it is appropriate to go suddenly from a scheme in which there is 
no charm parton distribution to one where the evolution and zeroth order coefficient function 
are such that the charm parton behaves much like any other parton, but that the zeroth order 
coefficient function should reflect the true physics near threshold. Indeed, there are several 
reasons why the above prescription is unsatisfactory. Usually at leading order in the expressions 
for structure functions involving only light quarks we include just the zeroth order coefficient 
function, and thus the expression is the summation of a s (Q 2 )ln(Q 2 ) terms only. Combining 
both Cf^ VF (Q 2 /m 2 ) and VF (z, Q 2 /m 2 ) with lowest order parton distributions we see that 
only the former really contributes to the LO expression, the first being down by a power of 
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\n(Q 2 ) at large Q 2 . However, keeping the simple form of VF above leads to an abrupt rise 
of the charm structure function as soon as Q 2 > m 2 , which is totally incompatible with the 



smooth behaviour in W 2 required physically. Thus, previous implementations |TIJ of a VFNS 
have included both both Cf) VF (Q 2 /m 2 c ) and C« VF (z,Q 2 /m 2 c ) in their LO definition of the 
heavy quark structure function, and improved the smoothness by also changing the scale from 
/i 2 = Q 2 to something like /i 2 = m 2 + 0.5 Q 2 (1 — m 2 /Q 2 ) 2 . Thus, we seem to have the choice 
of either a complete lack of smoothness in the charm structure function near threshold or an 
unconventional, and strictly incorrect ordering of the expression. We also have the disquieting 
feature that none of the heavy quark CF's contain the correct threshold behaviour in W 2 for 
F£. The correct threshold behavior can only be obtained overall by a cancellation of sometimes 
large incorrect contributions between different terms. All this just emphasises that the above 
procedure is not sufficient to yield a unique expression for the structure function in the VFNS. 

In order to better reflect the true physics in our choice of coefficient functions let us now 
consider the evolution of F% . From eq(|l|) the LO expression in the FFNS for the In Q 2 derivative 
is simply 

^(x,Q 2 ,m c 2 ) a. dCp FF (Q 2 /m 2 ) 

dlnQ 2 = Y J^Q 2 ® 9n '=* {Q } ' (14) 

The corresponding expression obtained by differentiating the LO expression in the VFNS, for 
Q 2 just above m 2 , is 



dF^{x,Q 2 ,m 2 ) dCP VF {Q 2 /m 2 c ) 
dlnQ 2 ~ d\nQ 2 



c + (Q 2 



+ g Cf VF (Q 2 /ml) ® (Pi ® g nf=4 (Q 2 ) + P° ® c + (Q 2 ))- (15) 

Throughout we are dropping the rif label on a s to simplify the notation. The question of how 
to treat the change in a s from rif = 3 to = 4 is discussed in our long paper ||. Now 



for Q 2 = m 2 , the terms in eq(r5) involving c + vanish because of eq(^I|) and so demanding 



continuity of the evolution across the transition point immediately leads, from eqs(T3,I3), to 

VF (Q 2 /ml) ® Pi = dC ° 1)F ^ /ml) . (16) 

Generalising this relation to be the definition of C^ -* VF (z, Q 2 /m 2 ) at allQ 2 offers a series of real 
advantages. Firstly, from the phenomenological point of view, it guarantees a smooth passage 
for the charm structure function from Q 2 < m 2 to Q 2 > m 2 , by definition. Explicitly taking 
the derivative of eq(0) is straightforward, giving 



dC^ 



FF, 



dlnQ 2 



z(l-2z 2 ) _ „ s l 



{P ° giz)+ 2e-^^-me 2 z 2 ) v 

\(W 2 -Am 2 c ) (17) 



+ {-Atz(l - 32) + 16eV) ln^^^^) + (4ez(l - z))v 
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and it is easy to see that in the limit Q 2 — > oo, 

Wr - W (18) 

Hence, from eg (P^6|) , we see that VF (z, e) must indeed tend to the simple form z 5(1 — z) in 
this limit, consistent with eq(|). Also, since FF (z, e) contains the factor 9{W 2 — Am 2 ) so 
does its In Q 2 derivative, thus ensuring the correct threshold behaviour in W 2 for CjP^ VF and 
in turn for F| at LO. Futhermore eq(^) allows the gluonic CF in the VFNS to be written as 



C^ F (z,Q 2 / m t) = C^ FF (z,Q 2 / m 2 c ) - ^^n*^ ln (S) ; 



(19) 



I) FF 

Finally, from eqs(|3| JT^JT8|) we see that as Q 2 — > oo, VF \z,Q 2 /m 2 ) does indeed tend to 



and so now VF has the same threshold behaviour as C^> FF 



the correct asymptotic MS limit given by eq(0). 

The extension of this procedure to any arbitrary order, i.e. continuity of the derivative 
in the gluon sector, is described in full in [Q. (Using our general approach continuity of the 
derivative in the singlet quark sector is not possible, but this has negligible effect in practice.) 

The only price that may seem to have been paid in exchange for obtaining a method which 
is far closer to true theoretical consistency and a much improved phenomenological description 
is the task of solving eq([TB|) for the LO charm CF, Cf^ VF (z, Q 2 /m 2 ). Note however that rather 
than CjP' VF appearing explicitly in the second term of eq(|l^), that term is now replaced by the 
simpler term in eq (|19"D and so we are just left with the term where CjP^ VF is convoluted with 
the charm density c + . Calculating this expression turns out to be easier than than deriving the 
form of CjP^ VF itself, the convolution being given by 



C? )FF (e)®c + (Q 2 ) = -Tdz 

J X 



dCp FF (z,e)fx\ 2 dc + (x/z,Q 



d\nQ 2 \zj d(x/z 

o r dC^ FF (z,e) x 
3/ dz « dlnQa -c + (x/z,Q) 



dz " V / dz'r(z')^c + (x/zz',Q 2 ) (20) 
dmQ A Jx/z zz 

where xq = (1 + 4e) _1 and r(z) is given by 

/ \ iT 1\ 3 . /v7 , 1\" 

r(z) = z 2 cos( — m-JH — ^ sm ( — hi -J . 

Thus we can calculate the charm structure function at LO, taking for Q 2 < m 2 



FF I 



(21) 



F^^ FF (x,Q 2 ,m 2 c ) = ^C^ FF (Q 2 /m 2 c ) ® g nf=3 (Q 2 ), (22) 
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and for Q 2 > m 2 

F c 2 ^ VF {x,Q 2 )=F^ FF {x,m 2 c ) + C^ VF (Q 2 /m 2 c ) c + (Q 2 ), (23) 

where the constant term becomes almost negligible for Q 2 > Am 2 . In the above, the partons 
are evolved only via the LO DGLAP equations. 

For the NLO case, the partons are evolved via the full NLO equations. The calculation of 
F^ FF for Q 2 < m 2 adds in the 0(a 2 ) contributions evaluated by Riemersma, Smith and van 
Neerven |T2[[] and in our VFNS we include the contributions from VF (Q 2 /m 2 c ) ® g nf =A and 
the extremely small contribution from C^p VF (Q 2 /m 2 ) ® £ n/ =4> where the coefficient function 
in the latter is determined by continuity of the evolution in the gluon sector at NLO. 

We also briefly discuss the treatment of the longitudinal structure function F£ since it 
is an important example. The 5-function CF for F 2 C of eq([H|), which previous procedures 
used, contains a contribution from F£ - corresponding to the 4m 2 /Q 2 term. Thus in those 
descriptions, _F£ in the VFNS contains a zero-order part, leading to a mis-match with the FFNS 
which starts at order a s with no \n(Q 2 ) enhancement. This implies rather odd Q 2 behaviour 
near threshold, see [O]. In our approach we do not require this zeroth-order contribution to 
F[ and so this problem is avoided. Our prescription ensures a smooth transition as in the F% 
case, see [[|. 

4 Results 

The results we present use a set of partons which are obtained from an analysis of DIS data 
in which charm and bottom flavours are treated by the new procedure we have discussed here. 
Only data on the full F 2 structure function are included in the fit, not the data on F% though 
we shall compare with those data below. The fits obtained are slightly sensitive to the value of 
m c and we find m c ~ 1.35 GeV gives the best fit. The value of m& was taken to be 4.3 GeV. 
We have evidence already that our more theoretically consistent treatment of heavy flavours 
improves the overall quality of the fits to F 2 ||. We plan to incorporate this procedure into a 
new (MRS type) global analysis where new data from hadronic collisions will be included and 
together provide updated parton distributions. 

In fig. 1 we show the result of the LO calculation of F% for two values of x. This demonstrates 
the success of our key result, namely the new form of the charm quark CF, which ensures a 
smooth transition from the description at low Q 2 in terms of the FFNS to high Q 2 in terms of 
the ZM-VFNS. In contrast to previous approaches JTU] this has been achieved without resorting 
to an unusual definition of LO structure functions (i.e. one in which the LO ZM-VFNS limit is 
not actually reached because 0(a s ) coefficients functions are included) or to any cancellation 
between large individual contributions. 

Fig.2 shows the similar situation for the NLO case compared with the NLO ZM-VFNS 
result and the continuation of the NLO FFNS expression. In the NLO case, the new description 

1 We are grateful to them for providing the program to compute these corrections 
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remains close to the continuation of the FFNS result to larger Q 2 than in the LO case. This 
is what is expected since, at the higher order, the FFNS includes more ln(Q 2 /m 2 ) terms which 
are resummed in the VFNS. 

In fig. 3 we show the comparison of the charm structure function F£(x, Q 2 ) resulting from 
our NLO analysis with all available data. The data at intermediate x values comes from EMC 



13] measurements of inclusive muons, the new data from HERA is by the HI U and ZEUS 
@ collaborations measuring D and D* cross sections. We show curves for three values of the 
charm mass and, as expected, the sensitivity to m c is strongest at low Q 2 . The values of F 2 C 
there from the EMC measurements suggest the higher value m c = 1.5 GeV is favoured. Clearly, 
precise measurements of F% in this region would provide very accurate estimates of the charm 
mass. 

5 Summary 

In this paper we draw attention to a new formulation for describing heavy flavour structure 
functions in DIS. While the concept of the various schemes discussed in section 2 is not new, 
the realisation (discussed in section 3) that beyond all orders consistency and order-by-order 
matching of the charm structure function across the transition from the FFNS to the VFNS 
extra constraints are required to uniquely determine all the relevant CF's is new. In fact, previ- 
ous attempts to describe charm evolution suffered from both theoretical and phenomenological 
deficiencies. By demanding that the order-by-order evolution of the charm structure function 
should be continuous as well as _F 2 C itself, we improve both the theoretical and phenomenological 
description of charm (and bottom) production in DIS. This successful procedure, based on an 
examination of the structure function expansion in a systematic order-by-order way, means that 
the procedure has a natural generalisation to any given order and also to any given physical 
process. 

We have demonstrated that the observed charm structure function is very well described 
over a wide range of x and Q 2 including especially the threshold region and that, consistent with 
the relatively large importance of charm at low values of x, fits to the total F2 are improved. 
For full details of our analysis for the general situation of heavy flavour production in DIS we 
refer the reader to ref|J. 
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Figure captions 



Fig. 1 Charm quark structure function, -F 2iC (x, Q 2 ) for x = 0.05 and x = 0.005 calculated using 
our LO prescription, our input parton distributions evolved at LO and renormalization 
scale fi 2 = Q 2 . Also shown are the continuation of the LO FFNS expression and the ZM- 
VFNS expression both calculated using the same parton distributions and same choice of 
scale. 

Fig. 2 Same as fig. 1 but with NLO prescriptions and NLO parton distributions. 

Fig. 3 Our prediction for F£(x, Q 2 ) using our NLO prescription, the NLO partons obtained from 
our global fit and three different values of m c compared with data from EMC |T3|], HI ]]J, 



ZEUS H] (from the 1994 run) and preliminary data from ZEUS (1995 run) WA . 
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